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Abstract 

The Telescope Array's Middle Drum fluorescence detector was instrumented 
with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. 
The data observed by Middle Drum in monocular mode was analyzed via the 
HiRes-1 profile-constrained geometry reconstruction technique and utilized the 
same calibration techniques enabling a direct comparison of the energy spectra 
and energy scales between the two experiments. The spectrum measured using 
the Middle Drum telescopes is based on a three- year exposure collected between 
December 16, 2007 and December 16, 2010. The calculated difference between 
the spectrum of the Middle Drum observations and the published spectrum 
obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy 
scale to be transferred to Middle Drum. The HiRes energy scale is applied 
to the entire Telescope Array by making a comparison between Middle Drum 
monocular events and hybrid events that triggered both Middle Drum and the 
Telescope Array's scintillator Ground Array. 

Keywords: UHECR, cosmic ray. Telescope Array, energy spectrum. High 
Resolution Fly's Eye, monocular, hybrid, HiRes 



1. Telescope Array 

The Telescope Array (TA) is the largest cosmic ray experiment in the north- 
ern hemisphere. It was designed to help resolve physics differences between the 
High Resolution Fly's Eye (HiRes) experiment, the Akeno Giant Air Shower 
Array (AGASA) 0, and the Pierre Auger Observatory [3] . TA consists of three 
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Figure 1: Map showing the location of Telescope Array relative to Salt Lake City and Dugway, 
Utah (the location of the High Resolution Fly's Eye). The route from Salt Lake City to Delta 
is 136 miles. 

HiRes-like fluorescence telescope stations overlooking 507 AGASA-like scintil- 
lator surface detectors (see Figures [1] and [2]) . The Surface Detector (SD) array 



was deployed in a square grid with a 1.2 km separation, covering ~ 700 km^ [l| 



Each SD unit is composed of two layers of 3 m^ x 1.2 cm scintillating plastic 
sheets separated by a thin steel sheet. The light from each layer is collected by 
wavelength-shifting optical fibers and directed into separate Photo-Multiplier 
Tubes (PMTs). 

Three telescope stations view the sky over the scintillator array. The north- 
ernmost fluorescence station, known as the Middle Drum site, consists of 14 
telescopes refurbished from the HiRes experiment's HiRes-1 site. These were 
deployed between November, 2006 and October, 2007 and were arranged to view 
~ 120° in azimuthal and 3° — 31° in elevation. Compared to HiRes-1 the 
Middle Drum site has only 1/3 of the azimuthal coverage but observes twice 
the elevation, as it was deployed into two rings, each covering 14° in elevation. 
Each telescope unit uses sample-and-hold electronics with a 5.6 fis gate. Each 
telescope camera consists of 256 PMTs covered with an ultra-violet band-pass 
fllter. Descriptions of the Black Rock and Long Ridge telescope stations were 
described by Tokuno [3]. 

The goals of the Middle Drum spectral analysis are three-fold. The primary 
goal of this analysis is to determine the flux of particles using the same cali- 
bration and analysis processing tools used to produce the monocular spectrum 
from the HiRes-1 data. The second goal is to compare the spectra measured 
by the Middle Drum detector with that of HiRes-1. Since the telescope units 
used in both of these detectors are composed of the same equipment, the re- 
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Figure 2: Map of the Telescope Array detectors. The telescope stations (Middle Drum, Long 
Ridge, and Black Rock Mesa) are indicated by blue triangles. The scintillator detectors are 
indicated by black dots and the Central Laser Facility by a red square. The ground array 
occupies about 700 km^ west of Delta, Utah. 



4 



suits of this comparison produce a direct link in the energy scale between these 
two experiments. Finally, by comparing events observed by Middle Drum and 
any of the other TA detectors, the energy scale of the entire Telescope Array 
experiment can be compared to that of the HiRes experiment. In this paper, 
this comparison is performed between the geometries of the events observed by 
Middle Drum and reconstructed using the monocular technique to those events 
that triggered both Middle Drum and the SD array and analyzed using a hybrid 
technique. 

2. Event Reconstruction and Selection 

The Middle Drum data and Monte Carlo events (described in section [3]) 
were processed using the same programs created for HiRes- 1 analysis Q. The 
only changes made were for the location and pointing directions of the tele- 
scopes. The HiRes-1 analysis was unique in that there was limited elevation 
coverage and a traditional monocular reconstruction could not be performed on 
the data. Instead, a combined geometrical-profile reconstruction was developed 
by Abu-Zayyad |7| which increased the resolution of the observed showers. This 
technique was not required in the analysis of Middle Drum data since the detec- 
tor observed longer-track events due to the increase in elevation angle coverage, 
however, it was used for consistency. 

As at HiRes, lasers are used for light-attenuation calibration, aerosol mea- 
surements, and relative-timing variances between the three fluorescence detector 
sites. Most of the events observed by the Middle Drum detector belong to these 
calibration lasers which are primarily removed by only processing those events 
that are downward-tending, since the lasers are fired in either upward or hori- 
zontal directions. Some of these laser shots appear to be downward-going events 
due to preliminary calculations using the timing and pointing directions of the 
triggered tubes. These are removed using the GPS trigger time-stamp and the 
GPS measured site positions. 

After filtering out laser events, most of those events that remain are due 
to electronic noise triggers, airplanes, and niuons that pass through the cam- 
era's PMTs. These are removed by determining a correlation between the time 
and geometrical pattern of the triggered tubes. Triggered tubes are clustered 
into groups of three or more tubes with difference limits on the trigger-time of 
2.0 /is and the pointing-direction of 1.2° from the previously triggered tube. 
These clusters are then combined into a single event-track from which a shower- 
detector plane (SDP) is determined. The tubes in a track are then iteratively 
checked and removed if greater than 3 RMS deviations away in either time or 
angle from the mean Q. 

The Middle Drum data and Monte Carlo simulations are reconstructed in 
monocular mode with the geometry determined by the equation 




(1) 



5 



where ti and Xi a-re the respective trigger time and pointing direction of tube 
i, Rp is the impact parameter of the shower with respect to the detector, if) is 
the angle the axis of the shower makes with respect to the direction of the core 
impact position around the detector, and to is the time the shower is calculated 
to be at Rp. 

Theprofile of the shower is calculated using the Gaisser-Hillas parameteri- 
zation [8| 

X- Xq 



Ne{x) 



Xo 



(2) 



where Np{x) is the number of charged particles (measured from the signal 
strength) at a given slant depth, a:, in g/cm^; N^ax is the maximum number of 
secondary particles produced in the extensive air shower, located at X^ax] Xq 
is a fit parameter associated with the depth of the first interaction; and A is a 
fit paramater defining the width of the shower profile. 

To reconstruct the Middle Drum data, as was done for HiRes-1, the time- fit 
was constrained by the shower profile reconstruction. This was performed by 
setting A to a constant 60 g/cm^ and Xq to a constant —100 g/cm^ in order to 
constrain the width and initial depth of the shower. These constants are in good 
agreement with average simulated shower measurements [6] . An inverse-Monte 
Carlo reconstruction is then made by simulating showers that closely resemble 
the true event using the triggered tubes. This is performed by choosing a series 
of Xmax values for individual Monte Carlo events over all energies in the shower 
library. To determine a best-fit profile reconstruction, a comparison is made 
between the light signal actually observed to the one simulated for each tube 
considered in the reconstruction. This is effective since both the timing and the 
profile fits are only dependent on the trigger time and pointing directions of the 
tubes used in the reconstruction, which determine the slant depth of the shower 
that each tube is observing along the axis of the shower. 

Separate chi-square minimizations are then performed on the timing and the 
profile reconstructions for each of the constant X^ax values chosen. The timing 
chi-square is calculated by 



Xtime — ^2 ^« 



to + I — I tan ' 



(3) 



with the error, ai, determined by the time to cross the face of a PMT. The 
profile chi-square is calculated by 



Xprofile — ^ ( "l2 ) i^i ^i)"^ (4) 



where, as in the timing fit, the sum is performed over the tubes within 3 RMS 
deviations away from the shower-detector plane. The observed signal, 5*°, is also 
used to calculated the uncertainty, erf, which is estimated to be S° + S. The 
constant, S = 200, is obtained through adding in quadrature the sky noise and 
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electronic fluctuations. Details of the reconstruction codes can be found in the 
dissertation by Abu-Zayyad [7]. 

The optimal reconstruction is then determined by calculating a best com- 
bined chi-squared for each X^ax fit using 

(5) 

where x is the normahzed chi-square value calculated as 

X = Xfit X (6) 

Xmin 

where Xfit is the chi-square of from the fit for each X„iax and NDF„iin and 
Xmin are the number of degrees of freedom and the chi-square, respectively, for 
the smallest chi-square reconstruction. As mentioned previously, this innovative 
technique was developed to reconstruct HiRes-1 data which had a limited ele- 
vation coverage. Future analyses of Middle Drum data will include traditional 
monocular reconstruction techniques. However, this method was used in this 
current analysis to provide a direct comparison to the spectrum observed by 
HiRes-1. Additionally, this technique results in a better resolution and aper- 
ture than an unconstrained time fit, even for the longer tracks observed at the 
Middle Drum site. 

After the selection of candidate events is obtained, quality cuts are performed 
on the fully reconstructed showers to remove any event that exhibits anomalous 
behavior. These cuts were optimized for the the short shower tracks observed 
by HiRes-1 0. The HiRes-1 analysis was ideal for cosmic rays with energy 
greater than lO^^eT^ since they could be observed from farther away and would 
appear as short tracks in the lower elevations. These cuts are applied to Middle 
Drum events since the higher-energy events would still have shorter tracks and 
the overlapping energy range between HiRes-1 and Middle Drum could then be 
directly compared. Additionally, this gives a baseline to future analyses. Events 
are retained if: 

1. the event reconstruct well, as determined by 

• not rejecting too many off-plane tubes, 

• there are enough slant-depth bins to fit a profile, 

• a minimum is attained, and 

• the modified geometry still parameterizes the timing fit; 

2. the angular tracks are > 7.9°, so that there are enough triggered tubes to 
provide a reliable reconstruction; 

3. the shower depth into the atmosphere observed by the first tube used in 
the reconstruction is < 1000 g/cm^, so the fit is not focusing on the tail 
of the shower; 

4. the in-plane angle, ip, is < 120°, to make sure the detector is not over- 
whelmed with Cerenkov radiation; and 

5. the area of the mirror observing the fitted track (away from the mir- 
ror/tube edges) is > 0.9 m^, to ensure there is not a bias in the recon- 
structed signal strength. 
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3. Monte Carlo Simulation 



The energy-dependent aperture of the detector is the product of the effective 
area and the sohd angle of acceptance. This is calculated using the equations 



(yll7)o (master) = 271^ (Rj 



p max 



R; 



p mm 



) X (1 - COS^max) 



and 



An{E) = {An)Q 



Nr 



recon {^^MC— recon ) 



(7) 



(8) 



* gen {EMC-gen) 

where Rp is the distance of closest approach of the shower, 9 is the zenith 
angle of the shower, Nrecon is the number of events reconstructed with energy. 



E 



MC —recon ; 



and Ngen is the number of events generated with energy Emc- 



Counting the number of events reconstructed at an energy folds in the detector 
bias into the energy spectrum calculation. Alternatively, the detector efficiency 
at a certain energy can be calculated by replacing Nrecon(,EMC -recon) with the 
number of events retained with a certain generated energy, Nrecon{EMC-gen)- 
The aperture of the Middle Drum detector has been calculated to be ~ 60% 
that of HiRes-1 for reconstructed energies of 10^^ ° eV, falling linearly to 40% 
at lO^o o eV (see Figure O. 
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Figure 3: The calculated Middle Drum aperture compared to that for the HiRes-1 detector. 

The CORSIKA-simulated shower library used by Middle Drum was the same 
generated for HiRes, using QGSJETOl as the hadron interaction model 
These showers were thrown with an isotropic distribution and consisted of ^ 10 x 
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the exposure of the Middle Drum cohected data. The Monte Carlo simulated 
only proton events between 10^^-^ eV and 10^^ ° eV using values as measured by 
HiRes below the GZK cutoff [l^. A spectral index of 3.25 was used below 
■|^q18.65 gY g^jjjj 2.81 above. The spectral set was thrown without simulating the 
GZK suppression ll| 12|. The lower energy range was thrown out to a range 
of 25 km from the telescope site, well beyond where the detector is incapable of 
triggering on the fluorescence light of a 10^*-^^ eV cosmic ray shower. The higher 
energy range was thrown out to 50 km. The simulated showers of both energy 
regions were thrown with a maximum zenith angle of 80°. The CORSIKA 
output is fed into the detector Monte Caro resulting in events which look exactly 
like real data and are subjected to the same reconstruction programs and quality 
cuts. 



3.1. Data-Monte Carlo Comparison 

To verify the adequacy of the Monte Carlo used for the aperture calculation 
and to confirm that the Monte Carlo closely models the real data, it is important 
to compare the distributions of the reconstructed Monte Carlo events and the 
data. These are shown in three energy ranges (io^^ "~^®'^ eV, loi®-^"!^ ^ eV, 
and > 10^^ " eV) in order to demonstrate that the Monte Carlo has the same 
energy evolution as the data. The variables chosen for this comparison are those 
that directly determine the aperture: the impact parameter, Rp (Figure SJ; the 
shower zenith angle, 9 (Figure [5|); and the shower azimuthal angle, 4> (Figure 

El). 

The impact parameter distribution directly determines the effective area of 
the aperture, and, as expected, the mean of the distribution increases along 
with the spread, or RMS, as the energy increases. The zenith and azimuthal 
angles directly determine the solid angle of acceptance. For all three parameters, 
the (black) data points and (red) Monte Carlo histogram distributions are in 
excellent agreement. 

It should also be noted that since the Middle Drum telescopes are pointing in 
the South-East direction, and that there is a quality cut removing many of those 
events that are pointing towards the detector, there is a depletion observed in the 
azimuthal distribution in this direction. This variance decreases with increasing 
energy since the impact parameter moves farther away from the detector and, 
therefore, there are fewer showers pointing above the 120° limitation. 



3.2. Resolution 

Resolution plots indicate how well the detector simulation and reconstruc- 
tion programs perform by comparing reconstructed values to generated values 
in Monte Carlo simulated events. The three primary parameters that show 
the quality of the reconstruction are the impact parameter {Rp) and the in- 
plane angle {ip) obtained from the geometrical reconstruction, and the energy, 
obtained from the profile reconstruction. These are determined for the same 
three energy ranges as the data-Monte Carlo comparisons to show trends in the 
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Figure 4: Comparison of the data and Monte Carlo distributions of the impact parameter, 
Rp, in three energy ranges: IQIS O-IS-S j^giS-S-lS O eV, and > lO^'' ^ eV. The Monte 
Carlo (red histogram) is in excellent agreement with the data (black points with error bars). 
The number of entries indicates the number of data events observed in that energy range. 



10 



18.0 <= log^gER3,„„<18 


.5 




1 ' ' ' ' 


1 ' ' ' ' 


1 ' ' ' 


TAMD 












Entries 3000 
Mean 36.57 
RMS 14.34 


250 

200 

^ 150 

100 

50 
n 














T 

i 


T 

l-BH 


T T 


: m 1 1 






















m-< I--- 
1 


1 

T 

i-u-L^ 










T i-i 

...l-jH.-^.t.. 








' ■ 








_ H 


u 










-i 








.... 


.... 


.... 


.... 


.... 




■ ■ ■ 1 


. . . .- 



10 20 30 40 50 60 70 80 90 




19.0<=loa E <20 


.5 






1 ' ' ' ' I ' ' ' ' 


|— I— I— I 


TAMD 




^.|0 necon 










Entries 191 
Mean 40.9 
RMS mnfi 


7U 

18 

16 

14 

z 12 

10 

8 

6 

4 

2 
n 




























































H 
























H 


H 






Hi 


i-i 1 












H 


H 






_Fi 


1- 


■-1 


■H 


H 














=- T--H 


i-t 


H 






1 




H 






















= HB-1- 












H 


!-■-! 

I4|-H 


I-I" 








^ 1^ 

■i-B-fS?-- 




1-Bh 














TLL-T 




















1- 










.... 


.... 


.... 


1 I 1 1_ 

.... 1 .... 1 ... In 


T _ 



en 



Figure 5: Comparison of the data and Monte Carlo distributions of the zenith angle, 9, in 
three energy ranges. The Monte Carlo (red histogram) is in excellent agreement with the data 
(black points with error bars). The number of entries indicates the number of data events 
observed in that energy range. 
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Figure 6: Comparison of the data and Monte Carlo distributions of the azimuthal angle, (p, 
in three energy ranges: iQl^ O-lS-S eV, 10l8-5-l9.0 g^^^j ^ 101^ ° eV. The Monte Carlo 
(red histogram) is in excellent agreement with the data (black points with error bars). The 
number of entries indicates the number of data events observed in that energy range. The 
moan of this distribution indicates that most events are pointing away from the detector. 
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reconstruction. With increasing reconstructed energy, the geometrical parame- 
ters show a trend of improving resolution (see Figures [7] and [8]) . For aU energy 
ranges, the energy resolution is on the order of 20% (see Figure 



4. The Energy Spectrum 

The measured energy flux spectrum includes data collected using the Middle 
Drum fluorescence telescope station between December 16, 2007 and December 
16, 2010 (see Figure ITO)) . The spectrum only includes data collected on clear, 
moonless nights with minimal cloud cover in the view of the detector for re- 
liable reconstruction. This amounts to ^ 2400 site-hours of data collection, 
corresponding to a 9% duty cycle. Multiplying this on-time with the aperture 
determines the Middle Drum exposure to be ^ 1/3 that of the final HiRes-1 
exposure. 

After three years of collecting data, 3859 events were observed. For each 
energy bin in which Middle Drum has observed events, the average number 
of events is ~ 32.5% that observed by HiRes-1. This is consistent with the 
Middle Drum exposure calculation. As was mentioned previously, an inverse- 
Monte Carlo technique is used in order to determine the energy of the shower. 
The Monte Carlo shower library, parameterized by the Gaisser-Hillas equation 
(see Equation [5]), is sampled for similar X^ax values and projected along the 
calculated geometry. The signal for each slant-depth bin of the simulated shower 
is then compared to the observed shower and chi-square values are calculated for 
a series of psi angles. The minimum chi-square value for a combined geometry 
and profile is determined to be the best-fit reconstruction. The showers are 
then distributed into tenth-decade energy bins from 10^* ° — 10^^ ° eV. The raw 
energy distribution for the data is shown in Figure 1111 

The flux is calculated by combining the number of events and the exposure 
per energy bin using the equation 



where ni^E) is the number of events in a given energy bin, E] AE is the width 
of the energy bin; Ail{E) is the energy-dependent aperture calculated from 
Equations [7] and E] and Aton is the on-time of the detector. This flux is often 
multiplied by the cube of the energy to flatten the spectrum in order to more 
clearly show the subtle features of the flux of these particles. Figure [1^ shows 
the spectrum as determined from the Middle Drum data overlaid with that from 
the two HiRes detectors' monocular reconstructions [oj. These are in excellent 
agreement in both normalization and shape. A consistency between the spectra 
measured by the Middle Drum detector and the HiRes-1 detector is determined 
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Figure 8: Resolution in the measurement of the in-plane angle, ^, of Monte Carlo simulated 
events shown in three energy ranges: iQis o-iS-S eV, iQiS-S-i^ o eV, and > lOi"" eV. The 
gaussian fit is used to determine the detector bias and resolution. 
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Figure 9: Resolution in the measurement of the energy of Monte Carlo simulated events shown 
in three energy ranges: iQis o-is.s gy, ioi8.5-i9.o gy, and > ° eV. The gaussian fit is 
used to determine the detector bias and resolution. 



16 



TAMD Monocular On-time 




Year.Day (year) 

Figure 10: Integrated data collection time as a function of elapsed time for the Middle Drum 
fluorescence telescope site. 



where Jtamd and JniBes-i are the measured flux aud <Jtamd and (THiRes-i 
are the statistical uncertainties observed by the Middle Drum and HiRes-1 de- 
tectors, respectively (see Figure [T31)- This calculation only included those energy 
bins with at least seven events observed by the Middle Drum telescopes within 
the HiRes-1 spectral range. This difference shows that the flux measured by 
Middle Drum is within O.Scr of HiRes-1. They are consistent with the same flux 
level. 

Further, the Middle Drum spectrum can be quantitatively compared to the 
HiRes-1 spectrum by determining the between the flux measurements on a 
bin-by-bin basis. The value is calculated by summing over the square of each 
A J given in Equation [lOl The resuh is a x^/N.D.F. = 12.20/10 for aU of the 
overlapping bins and /N .D.F. — 4.47/5 for bins > 10^^ ° eV. Figure [T^ shows 
that the is dominated by the difference in the measured flux in the lO^^-^^ 
eV energy bin. 

Since Middle Drum measured the spectrum with the same equipment and 
calibration techniques and obtained the same result as HiRes-1, the HiRes-1 
energy scale is thus transferred to Middle Drum. Had the energy scale changed, 
the rapidly falling spectrum would have shifted by twice that increment. 
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TA Middle Drum Events 
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Figure 11: The raw energy distribution of events observed by the Middle Drum detector over 
the first three years of collection. The events are shown distributed in tenth-decade energy 
bins between lO^* " eV and lO^o o eV 
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Figure 12: The energy spectra multiplied by E^. The spectrum as determined from the Middle 
Drum data is shown by the black boxes. The spectra of HiRes-1 (upward, green triangles) 
and HiRes-2 (downward, magenta triangles) are shown for comparison. The three spectra are 
in excellent agreement in both normalization and shape. 
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TAMD HiRes-1 Energy Scale Difference 
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Figure 13: The bin-by-bin energy flux measurement difference between Middle Drum and 
HiRes-1. Only energy bins between lo^^-^^ and lo^^-^^ eV were compared since these are the 
only bins in the Middle Drum data where there are more than seven observed events within 
the HiRes-1 energy range. 
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5. Middle Drum Hybrid Geometry Comparison 



The transfer of the energy scale from the HiRes-1 spectrum to the Middle 
Drum spectrum creates a direct link between the HiRes and Telescope Array 
experiments. The next step to completely bridge the two experiments is to 
determine the energy scale between those events observed by the Middle Drum 
detector to those that also triggered the ground array. This Middle Drum 
monocular-hybrid comparison will then transfer the energy scale of HiRes to 
the rest of Telescope Array in future studies. 

The hybrid analysis begins by improving the geometrical reconstruction. 
Time and pulse height information of the triggered ground array scintillator 
detectors (SDs) are used to improve the time-versus-angle fit 13]. This is per- 
formed by calculating the SD core using the modified Linsley shower-shape 14 1 
to obtain a lateral distribution function (LDF) which is then used to constrain 
the monocular time-versus-angle fit performed for the FD. 

After the improved shower geometry is determined, the profile fit is per- 
formed using the inverse-Monte Carlo technique presented in this paper, how- 
ever, in this fit, the geometry determined above was not adjusted to scan for a 
better profile fit. The hybrid data selection cuts use a combination of Middle 
Drum and SD information. Events are retained if: 

1. the profile fit reconstructs well, as is determined in the monocular recon- 
struction; 

2. the geometry fit has a ^ jNDF < 7; 

3. the zenith angle is < 56°, providing a well- reconstructed SD core impact 
location of simulated showers thrown up to 60°; 

4. the SD calculated core must be within 500 meters of the SD boundary, so 
there is no bias in the LDF reconstruction; 

5. the SD calculated core must be within 600 meters of the shower-detector 
plane, so the shower track remains consistent between the two detectors; 

6. the angular track length is > 7.9°, to provide a reliable profile fit; and 

7. Xmax is observed by Middle Drum, for reliable composition studies. 

Since the boundary of the ground array begins ~ 7 km from the Middle Drum 
site, most of the monocular events with energy less than 10^^ '* eV fall outside of 
the ground array (see Figures 14(a) through |l4(d)[ ). Above this energy, roughly 



half of the events observed monocularly have core positions within the boundary 
of the ground array (see Figures [15(a)] and p(b)| ). 



Figure [16] shows a comparison of the energy for events reconstructed from 
the Middle Drum data in monocular mode to the same events reconstructed 
in Middle Drum-hybrid mode. Only events retained in both the monocular 
and hybrid analyses were compared. For those events with ^ Emono x Ehybr > 
2^018.5 ^Yie monocular and hybrid energies are in good agreement (see Figure 
[TT]) . This provides a direct link between the events observed by Middle Drum 
in monocular mode to those events that also trigger the ground array. A direct 
comparison is thus made between HiRes and all of the Telescope Array detectors. 
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MD Monocular Core Positions: 18.0 - 18.1 eV lUID Monocular Core Positions: 18.1 - 18.2 eV 




(c) 1018-2 - 1018-3 eV (d) 1018-3 - 1018-4 eV 



Figure 14: The core positions of the Middle Drum events observed and reconstructed in 
monocular mode (indicated by the points) for tenth-decade energy bins between 1018-*^ and 
10l8-4 eV. The locations of the fluorescence stations (BRM, LR, and MD), as well as the 
Central Laser Facility (CLF), are indicated by stars. The perimeter of the scintillator Surface 
Detector (SD) array is indicated with lines. At these lowest energies, virtually all core locations 
are outside of the SD array. 



22 




^ West East -> (1-2 km) West East ^ (1 .2 km) 



(a) - 1018-5 eV (b) 1018-5 _ iqIS-G eV 

Figure 15: The core positions of the Middle Drum events observed and reconstructed in 
monocular mode (indicated by the points) for tenth-decade energy bins between IQl* * and 
-|^q18.6 gY Compared to those events shown in Figure [T4l as the event energy increases, an 
increasing percentage of cores are observed within the SD array. 

6. Conclusions 

The Telescope Array's Middle Drum observatory uses refurbished telescopes 
from the High Resolution Fly's Eye experiment. A spectral measurement was 
made using the first three years of the Middle Drum data collection. Both 
the data and simulated events were analyzed monocularly using the profile- 
constrained geometry reconstruction technique that was developed for the HiRes- 
1 data. The energy and geometrical resolutions of the Monte Carlo simulations 
show good agreement between what was generated and what was reconstructed 
and the data-Monte Carlo comparisons are in excellent agreement between sim- 
ulated and real extensive air showers. The calculated Middle Drum energy 
spectrum is shown to be in excellent agreement with the spectra produced by 
the HiRes-1 monocular analysis with the difference between them less than the 
energy resolution of the Middle Drum reconstruction. The HiRes energy scale 
can now be transferred to the entire Telescope Array for further comparisons 
now under way. 
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Energy Difference TAMD Mono.-Hybr. 




18.2 18.4 18.6 18.8 19 19.2 19.4 19.6 19.8 20 20.2 



'OQlO ^Monocular (©V) 

Figure 16: Comparison of event energy for cosmic ray events reconstructed in monocular mode 
using data collected from the Middle Drum fluorescence telescope station (abcissa) versus the 
energy when reconstructed in hybrid mode (ordinate) incorporating information from the 
scintillator surface array. The solid line drawn indicates where the two measurements would 
be equal. There is excellent agreement between the measurements. 
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Figure 17: Spread of difference in event energy for cosmic ray events observed by the Middle 
Drum fluorescence telescope station reconstructed in hybrid mode incorporating information 
from the scintillator surface array to that reconstructed in monocular mode. The gaussian 
mean shows the overall difference between these reconstruction techniques, with an uncertainty 
indicated by the standard deviation. There is excellent agreement between the measurements. 
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